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Herein we report the discovery and SAR of a novel series of SARS-CoV 3CLpro inhibitors identified 
through the NIH Molecular Libraries Probe Production Centers Network (MLPCN). In addition to 
ML188, ML300 represents the second probe declared for 3CLpro from this collaborative effort. The X- 
ray structure of SARS-CoV 3CLpro bound with a ML300 analog highlights a unique induced-fit reorgani¬ 
zation of the S2-S4 binding pockets leading to the first sub-micromolar noncovalent 3CLpro inhibitors 
retaining a single amide bond. 

© 2013 Elsevier Ltd. All rights reserved. 


Coronaviruses (CoV) are enveloped, large plus-strand RNA 
viruses associated with mild to severe respiratory symptoms, 
including the common cold and the Severe Acquired Respiratory 
Syndrome (SARS)-CoV. Identified as the etiological agent 
responsible for the global pandemic in 2003, SARS presents an 
atypical pneumonia that during the first major outbreak led to pro¬ 
gressive respiratory failure in over 8000 individuals and about 800 
deaths by July of that year. 4 With the cooperation of leading 
nations, a rigorous public healthcare campaign was fortunately 
successful in controlling this outbreak. However, a reemergence 
of the SARS-CoV is considered a potential pandemic risk and new 
strains of human coronavirus continue to be identified. Since 
2003, two additional human coronaviruses, NL63 and HKU1, have 
been identified in patients and the viruses have been characterized 
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and found to be significantly less lethal than SARS-CoV. Most re¬ 
cently in 2012, a new SARS-like virus, designated the Middle East 
respiratory syndrome coronavirus (MERS-CoV), has been identified 
in 144 patients so far, 54 of whom died. 8 There is now evidence for 
person-to-person transmission of MERS-CoV. Now, nearly a dec¬ 
ade later, the possibility of another SARS-like pandemic appears 
even more palpable based upon the lethality and properties of 
the newly identified MEV-HCoV strain. Effective vaccines and small 
molecule antiviral agents to prevent or treat SARS-like infections 
still do not exist, thus tailored antiviral therapies are urgently 
needed in order to treat potential future outbreaks of SARs and re¬ 
lated human coronaviruses. 

The SARS and MERS coronaviruses encode two proteases, a pa- 
pain-like protease (PLpro) and a 3-chymotrypsin-like protease 
(3CLpro), in their genome that are essential for viral replication. 
The viral polyprotein is cleaved at three unique sites by PLpro 
and 11 unique sites by 3CLpro. Initial reports of 3CLpro inhibitors 
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in the literature focused on peptidomimetics, often four to five res¬ 
idues in length, bearing a reactive ‘warhead' group, such as an alde¬ 
hyde, halo-methyl ketone, or Michael acceptor at the terminus 
with several demonstrating a covalent interaction with the active 
site Cys-145 residue. 10 1 Until recently, the majority of efforts to 
develop nonpeptidic 3CLpro inhibitors also relied on ‘warhead’ 
based design strategies (Fig. 1, 1—5) 1 21 and a number of these 
nonpeptidic inhibitors achieved sub-micromolar activity. In the 
case of pyridyl ester 4, this potent nanomolar mechanism-based 
enzyme inactivator led to cell based inhibition below 10 pM in 
SARS-CoV infected Veto E6 cells. Recently, we reported N-(tert- bu- 
tyl)-2-(N-arylamido)-2-(pyridin-3-yl) acetamide 6 (Fig. 1, ML188) 
and its X-ray complex with 3CLpro (PDB: 3V3M) as a rare example 
of a noncovalent SARS-CoV 3CLpro inhibitor of moderate molecular 
weight with good enzyme and antiviral inhibitory activity. ' Here¬ 
in, we describe the continuation of efforts to develop potent, non¬ 
covalent SARS-3CLpro inhibitors based upon a second chemical 
class of triazoles from our MLPCN screening campaign (7, Fig. 2) 
and progression of this lead series to a second generation probe 
ML300 (16e, Fig. 1) and beyond to arrive at sub-100 nM inhibitors. 
We propose from crystallography data that ML300 and related tri¬ 
azoles in this series inhibit 3CLpro via a novel mechanism of action 
and provide a new direction for additional noncovalent inhibitor 
design and refinement. 

Using a designed expression construct which produces the 
post-proteolytic and authentic 3CLpro dimer, a screen against the 
NIH molecular libraries sample collection (~293 K compounds) at 
the Scripps Research Institute Molecular Screening Center 
(SRIMSC) was undertaken. In addition to the diamide acetamide 
series represented by ML188 (6, Fig. 1), 2 ~ a related diamide series, 
represented by SID 24808289 (7, Fig. 2), was identified demon¬ 
strating a 3CLpro IC 50 of 6.2 pM and good selectivity versus PLpro 
(1C 50 > 60 pM) which is used as a control for cysteine-protease 
activity. Fortunately, quite early in the chemistry campaign an X- 
ray crystal structure of diamide 7 bound to 3CLpro was determined 
to 1.85 A resolution. A solvent accessible surface depiction of 7 in 
the 3CLpro active site along with a wall-eye stereo view with key 
contact resides and hydrogen bonding contacts in depicted in 
Figure 3. Interestingly, in contrast to the ML188-3CLpro crystal 
structure in which ML188 accommodates substrate sub-pockets 
in the enzyme active-site traditionally occupied by peptidomimet¬ 
ics, diamide 7 engenders an induced-fit complex resulting in a new 
surface dictated largely by a rearrangement of the Gln-189 and 



1 (2004) 17 
MAC-5576 
IC 50 = 0.5 pM 


2 (2005) 18 
Cinanserin 
IC50 = 5 pM 



3 (2006) 19 


K, = 11.1 fiM 



4 (2008) 20 
1C 50 = 30 nM 

SARS-CoV EC 50 = 6.9 pM 



5 (2008) 21 
IC 50 =13^M 

SARS-CoV EC 50 = 12.9 pM 
pdb code: 3V3M 


6 , ML188 (2013) 22 
IC 50 = 1.5 pM 


Figure 1. Representative nonpeptidic 3CLpro inhibitors utilizing a warhead and 
noncovalent mechanism of inhibition (1-6). 


P2-P4 P2-P1' 




7 (SID 24808289) 
3CLpro IC 50 = 6.2 pM 
PLpro > 60 pM 


16e (ML300) 

3CLpro IC 50 = 4.11 ± 0.24 pM 


Figure 2. Binding orientation and properties of MLPCN 3CLpro HTS hit 7 and 
evolved probe molecule ML300(16e). 




Figure 3. (A) Solvent accessible surface view of 7-3CLpro complex (PDB code: 
4MDS, PubChem SID 24808289); (B) X-ray crystal structure of 7 (capped sticks in 
orange carbon) with SARS 3CLpro in wall-eye stereo view with key residues and 
hydrogen bonds. 


Met-49 residue side-chains. 23 This induced fit accommodates the 
syn N-methyl pyrrole and anilido acetamide moieties of the 
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inhibitor within subpockets that can be characterized as S 2 -S 4 and 
S 2 —SF subpockets, respectively. 

Figure 2 schematically illustrates the inhibitor-active site inter¬ 
actions oriented in a manner similar as depicted in Figure 3. In 
addition to the P 2 -P 4 and P 2 -Pj' groups the inhibitor partially 
occupies the S 3 subpocket with a terminating 2-methylbutylamide. 
Key hydrogen bonding interactions can be found near the catalytic 
site with His-163 and the benzotriazole N-( 3) engaged in a key 
interaction, with an interatomic distance of 2.9 A. In addition a 
backbone Glu-166 NH interaction is evident with the central acet¬ 
amide oxygen (N-0 distance 2.8 A). 

Flexibility of the diamide scaffold (RotBon ~7) coupled with the 
observed induced-fit within the active site of 3CLpro presents an 
added challenge with respect to in silico inhibitor approaches. 
Thus, our structure-activity-relationship (SAR) studies focused ini¬ 
tially on three key areas within the diamide scaffold: (1) benzotri¬ 
azole replacements with alternative hydrogen bond acceptor 
functionality to interact with His-163, (2) acetamide modifications 
within the P 2 -P]' region, and (3) minimum pharmacophore dele¬ 
tion studies of the P 3 2-methylbutylmide. The P 2 -P 4 group was 
held constant for this investigation and based upon HTS and recon¬ 
firmation results (data not shown) the N-methyl pyrrole was re¬ 
placed with an equipotent 3-thienyl moiety. 

In parallel with efforts to obtain the 3CLpro-7 crystal structure, 
synthesis of first generation analogs to survey diversity of the ben¬ 
zotriazole unit were initiated using a modified version of our 4CC- 
Ugi strategy (Scheme 1) to allow for late stage azole introduction. 
Thus, Ugi reaction using f-butyl isocyanide, chloroacetic acid, 
thiophene-3-carbaldehyde, and N-(4-aminophenyl)acetamide pro¬ 
ceeded smoothly to give chloride 8, which could be isolated in 
good yield after chromatography. Displacement of chloride 8 with 
azole NH heterocycles provided 9a-c. Alternatively, displacement 
of 8 with sodium azide and subsequent Huisgen cycloaddition 
reaction with an appropriate acetylene furnished 1,2,3-triazoles 
lOa-c in good overall yield. 

Synthesis of P 2 -P]' amide analogs within the elaborated dia¬ 
mide were similarly prepared in an Ugi reaction using Boc-pro- 
tected 4-(amino) aniline (Scheme 2) as the amine component. 
Deprotection of 11 using trifluoroacetic gave aniline 12 which 
was reacted with a variety of carboxylic acid derivatives under 
HATU coupling conditions or reacted with an acid chloride or sufo- 
nyl chloride in the presence of TEA to give final examples 13a-l. 

Synthesis of P 3 truncated analogs began with reductive amina- 
tion using thiophene-3-carbaldehyde with either 4-bromoaniline 


NHAc 



NHAc 



9a-c 


NHAc 





Scheme 1. Synthesis of P! analogs 9a-c and lOa-c. Reagents and conditions: (a) 
MeOH, 50 °C, 4 h, 95%, (b) (i) NaH, HetNH, DMF, (ii) 9, DMF, 65-80%, (c) (i) NaN 3 , 
DMF, 100 °C pwave 30 min, 95%, (ii) acetylene (R = Ph, TMS), DCE, 120 °C 16 h, 85- 
98%, (iii) R = TMS, TBAF, HOAc, 0 °C-rt, 45%. Final library compounds were purified 
by UV prep or mass-directed prep HPLC. 


NHBoc 



12 13a-l 


Scheme 2. Synthesis of P 2 -Pi analogs 12 and 13a-l. Reagents and conditions: (a) 
MeOH, 50 °C, 4 h, 94%, (b) TFA, 95% (c) (i) HATU, DIPEA, DMF, RC0 2 H, 55-73% (d) 
RC0 2 C1 or RS0 2 C1, TEA, DCM, 51-64%. Final library compounds were purified by UV 
prep or mass-directed prep HPLC. 


R 



15a R= NHBoc 16a-k amide library 

15b R = Br 17a-e biaryl library 


Scheme 3. Synthesis of P 2 -PT analogs 16a-k and 17a-e within truncated series. 
Reagents and conditions: (a) NaHB(OAc) 3 , DCE, rt, 80% (b) benzotriazol-l-yl-acetic 
acid, HATU, TEA, DMF, rt, 74% (c) (i) 15a, TFA, DCM, 95%, (ii) HATU, DIPEA, DMF, 
RC0 2 H, 65-80%; RC0 2 C1 or RS0 2 C1, TEA, DCM, 90-95%; NaHB(OAc) 3 , RCHO, DCE, 
45-95% (d) 15b, Ar/HetB(OH) 2 , 1 M aq Na 2 C0 3 , 5mol% Pd(PPh 3 ) 4 , THF, 30-65%. 
Final library compounds were purified by UV prep or mass-directed prep HPLC. 


or Boc-protected 4-(amino) aniline to give intermediates 14a-b 
in good yield. Amide coupling with HATU using benzotriazol-l- 
yl-acetic acid installed the requisite P] groups to afford 15a-b. Ini¬ 
tial efforts focused on preparing the identical amide library 
prepared in the elaborated series 13 (Scheme 3). This was readily 
accomplished as before; Boc-deprotection of 15a followed by 
amide coupling or acylation/sulfonylation, gave 16a-k. Subsequent 
synthesis of a series of biaryls as amide replacements commenced 




9a 9b 

Inactive Inactive 


9c 

Inactive 


Cm 



10 a 10b 

Inactive Inactive 


10 c 

3CLpro IC50 = 11 pM 


Figure 4. Representative azole replacements (9a-c and lOa-c). 
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using a Suzuki cross-coupling with 15b and a variety of boronic 
acids to afford target molecules 17a-e. 

SAR for 1,3-azole Pt replacements (9a-c, Fig. 4) indicated a 
strict requirement for the 1,2,3-triazole unit; benzimidazoles 
9a-b and 2-methyl-l-imidozyl derivative 9c were uniformly inac¬ 
tive. Since the N-(3) nitrogen of 7 appeared to be involved in a 
hydrogen bond with His-163, it was somewhat surprising that 9a 
was not tolerated since the imidazole has the potential to maintain 
a N-(3)-His-163 hydrogen bond interaction. However, within the 
3CLpro-7 structure the catalytic Cys-145 residue is located within 
3.3 A of the N-(2) nitrogen, indicating potential for a weak 
hydrogen bond and/or dipole-dipole stabilization interaction. This 
potential interaction may perhaps be responsible for the 1,2,3-tri- 
azole preference. Interestingly, 4-phenyl 1,2,3-triazole 10c was tol¬ 
erated with an 1C 50 of 11 pM, suggesting additional avenues for 
optimization. Accommodation of the phenyl moiety of 10c within 
the active site St subpocket is not entirely clear at this time. Based 
on the 3CLpro-7 structure, Glu-166, Phe-140, and Glu-166 are pre¬ 
dicted to be within close proximity. Unsubstituted triazole 10a and 
trimethyl silyl triazole 10b were inactive, demonstrating the 
importance of maintaining a proper aromatic ring in this 
subpocket. 

Amide library 13a-l (Fig. 5) within the elaborated diamide ser¬ 
ies displayed a range of potency from moderate micromolar activ¬ 
ity (13a, 13b-d, 13f-g) comparable to the HTS hit 7, to weakly 
active or inactive. Cyclic and acyclic acetamide congeners related 
to HTS lead 7 showed consistent activities below 10 pM with 
branched i-propyl derivative 13d and cyclobutyl amide 13g having 
the greatest activity below 5 pM. Modification to sulfonamide 13b 
resulted in a three-fold loss in inhibition relative to acetamide 13a. 
The smaller cyclopropyl (13f) or larger cyclohexyl (13h) cyclic 
amide generally resulted in loss of inhibition. Incorporating a ste- 
rically hindered t-butyl amide 13e also led to a modest three-fold 
loss in activity. Lastly, aromatic and heteroaromatic amides (13i- 
k) in addition to iso-butyl carbamate 131 were weak or inactive 
as 3CLpro inhibitors. Collectively, these data appear to be consis¬ 
tent with the 3CLpro-7 structure whereby a short helix-loop-helix 
motif (Val-42-Ile-43-Cys-44-Thr-45-Ala-46) and a proximal p-turn 


NH 


2 





P2-P1' Exploration: 
amide, urea, sulfonamide 



12 

IC 50 > 20 nM 
poor curve fit 

O 


13a 13b 13c 

IC 50 = 7.72 HM IC 50 = 25.3 M-M IC 50 =6.9^iM 


13e 

IC 5 o= 22.5 jjM 




13f 

IC 50 =9.1 hM 


13d 

ICso=4.1 


13g 

IC5o = 3.8 jiM 


13h 

23% inh (100 |iM) 



Inactive 15% inh (100 (iM) IC 50 = 26 pM Inactive 


Figure 5. 3CLpro activity from library 13. 


element (Thr-24-Thr-25) define the edge of this pocket with min¬ 
imal volume for larger groups beyond acetamide 7. 

With limited success from the above S! and S 2 -Si' studies we 
turned to P 3 truncation to potentially identify a minimum pharma¬ 
cophore to reduce overall MW and improve ligand efficiency (LE). 24 
Examination of the P 3 group within the 7-3CLpro structure 
suggested this group was unfavorably solvent exposed relative to 
the t-butylamide-S 3 interaction found within the ML188-3CLpro 
structure. 22 Initial efforts led to 16a-k (Table 1). Gratifyingly 
truncated amides proved to have comparable activity versus the 
elaborated diamide counterparts (Fig. 5 see 13c-d, 13f-g, vs Table 1 
see 16a-c, 16e-f). Interestingly, truncated series 16 appeared to 
better tolerate larger substituents, perhaps suggesting additional 
changes in the shape of the active site within this subpocket. For 
example cyclohexyl amide 16g was found to be a weak inhibitor 
and similarly carbamate 16i had moderate inhibitory activity of 
10.3 pM while its related counterpart 131 (Fig. 5) was inactive. 

At this stage in the project with efforts focused on P 3 truncated 
analogs bearing a putative S 2 -S 1 ' interaction, we elected 16e for 
further characterization and probe declaration (ML300, Fig. 6). Rel¬ 
ative to probe ML188 (6) and the equipotent diamide 13d, ML300 
proved to offer improvements in several areas (Fig. 6). SARS 3CLpro 
inhibitor ML300 is ~100 amu lower MW (MW = 431) relative to 
13d with moderate ligand efficiency (LE). Deletion of the lipophilic 
P 3 group reduces cLogP an order of magnitude (cLogP = 3.2) and 
thus greatly improves ligand-efficiency-dependent lipophilicity 
(LELP) versus ML188 and 13d. Probe molecules ML188 and 
ML300 were evaluated in an in-house 2: in vitro DMPK panel 
including plasma protein binding, P450 enzyme inhibition, and 


Table 1 

3CLpro activity 16a-k, 17a-e 


16a-k 

17a-e 


Cmpd R IC 50 a Cmpd R IC 50 a 



a IC 50 are the average of three independent determinations and represent a 
co-efficient of variation (CV) < 0.10. 
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SID 99350510, MW=433 
3CLpro IC 50 = 1.5 pM 
LE=0.25 

cLogP=4.6, LELP=18 



13d (racemic) 

SID 93373783, MW=532 
3CLpro IC 50 = 4.1pM 
LE=0.21 

cLogP=4.0, LELP=19 



16e (ML300) 

SID 99289112, MW=431 
3CLpro IC 50 = 4.11 ± 0.24 pM 
LE=0.24 

cLogP=3.2, LELP=13 


In vitro DMPK and ancillary pharmacology: 
ML188 

PPB Fu (h, r): 3.3%, 2.6% 

CLhep (h, r): 20, 69 mL/min/kg 
CYPs (pM): 

>30 (1A2), 9.1 (2C9), >30 (2D6), 1.7 (3A4) 
Eurofins Profiler Screen: 

No significant activity 


ML300 

PPB Fu (h, r): 4.8%, 10.6% 

CLhep (h r): 20, 67 mL/min/kg 
CYPs (pM): 

7.7 (1A2), 7.2 (2C9), 8.4 (2D6), 4.6 (3A4) 
Eurofins Profiler Screen: 

Melatonin MT1 75% inhibition (10 pM) 


Figure 6. Profiles of 3CLpro inhibitors 6 (ML188), 13d, and 163 (ML300). 


intrinsic clearance using liver microsomes (Fig. 6). Both ML188 and 
ML300 possess good free fraction with ML300 being superior (1.5 
and 4.0-fold improved human and rat fraction unbound, respec¬ 
tively); however, intrinsic clearance (CL HE p normalized to liver 
blood flow, Qh = 21 mL/min/kg, Q, = 70 mL/min/kg) indicates 
ML188 and both ML300 are both predicted to be highly cleared. 
ML188 and ML300 possess modest P450 enzyme inhibition, with 
ML300 maintaining 5-10 pM activity across four major CYP 
enzymes (Fig. 6). Probe ML300 was found to be highly selective 
in a Eurofins lead-profiling screen 26 with only modest activity for 
melatonin MT, receptor in a radioligand binding assay. 

Based on the absence of key hydrogen bonding interactions of 
the P 2 -P 1 ' amide of 7 with the 3CLpro active site, in addition to 
the poor metabolic instability and CYP profile of ML300, we opted 
to explore more diverse amide replacements as a means to im¬ 
prove metabolic stability, P450 activity, and 3CLpro inhibitory 
potency. Initial efforts identified representative N-methyl (16j) 
and N-benzyl (16k) anilines with potency comparable to probe 
ML300 (16e). The lack of activity for benzamide 16h versus the 
reduced benzylamine 16k is striking and indicates the enhanced 
flexibility of the N-benzyl group is permitting a productive interac¬ 
tion to occur where previously aromatic amides were not tolerated 
(see Fig. 5, 13i-j). A subsequent survey of biaryls was explored 
(Table 1, 17a-e) and on the basis of the 3CLpro-7 X-ray, we 
targeted 3-pyridyl (17b-c) and 4-pyrimidyl (17e) heterocycles as 
means to potentially engage a side-chain interaction from the 
hydroxyl groups of Thr-24 or Thr-25. These modifications afforded 
inhibitors with micromolar activity and in the case of 2-methoxy- 
pyridyl 17c submicromolar activity (IC 50 = 700 nM). Unexpectedly 
the parent simple phenyl biaryl 17a proved to have a major impact 
on activity with a ~7-10-fold increase relative to 17b-c. 3CLpro 
inhibitor 17a represented the first sub-100 nM inhibitor for the 
series and to our knowledge one of the most potent nonwarhead 
based SARS 3CLpro inhibitors to date. At this time inhibitor 17a 
is relatively unoptimized and thus current efforts are focused on 
targeted biaryl congeners to understand DMPK, cellular activity, 
as well as potential broad spectrum activity against other corona- 
virus strains including MERS-CoV. 

In summary, we have described the identification and binding 
orientation and interactions for a second class of diamide SARS 
3CLpro inhibitors, culminating in probe ML300 and subsequently 
improved inhibitors such as 17a, which possess LE > 0.3 and an 
LELP approaching 10. The X-ray crystal structure of HTS hit 7 
bound to 3CLpro was instrumental in guiding optimization and 


the induced-fit of this inhibitor 3CLpro complex illustrates the 
challenges of divergent SAR and the limitations of virtual based 
screens. The four component Ugi reaction was utilized once more 
to rapidly generate SAR for the putative P 2 -Pi' and Pt subgroups. 
Importantly, P 3 truncation was possible for this triazole series of 
3CLpro inhibitors, allowing for significant MW reduction without 
diminishing potency. Collaborative efforts in these laboratories 
continue towards the identification active inhibitors within the 
truncated biaryl class. Integrated efforts with DMPK assessment 
continue in order to improve intrinsic clearance and diminish 
P450 activity, which are issues to be addressed within the series 
prior to in vivo proof-of-mechanism studies. ML300 is an MLPCN 
probe and is freely available upon request. 
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